Spectra of graphs and semi-conducting 

polymers 

00 ' 

f^ ' Philipp Schapotschnikow and Sven Gnutzmann 

o: 

[•-[- . Abstract. We study the band gap in some semi-conducting poly- 
mers with two models: Hiickcl molecular orbital theory and the so- 

^J . called free electron model. The two models are directly related to 

,——1 ' spectral theory on combinatorial and metric graphs. Our numeri- 

rS^ . cal results reproduce qualitatively experimental results and results 

r^' from much more complex density-functional calculations. We show 

r^ , that several trends can be predicted from simple graph models as 

"^ ' the size of the object tends to infinity. 

g. 

1. Introduction 

K^ I An important pattern in organic chemistry are delocalized electron 

^ I systems, also called conjugated 7r-electron systems. Most dyes owe 

^ ■ their color to such structures. When the length of a conjugated elec- 

^ ! tron system becomes large, electric conductivity and luminescence are 

'""' I possible. 

f^^ ' This is achieved in practice by chemically connecting a fragment with 

O ■ delocalized electrons to its copy, and repeating this procedure many 

^ ! times. The original fragment is called monomer; if the number of copies 

'k>( I is definite and rather small (up to few tenths), the resulting product 

5^ ■ is called oligomer; if the number of segments is large and indefinite, 

the resulting product is called polymer. Semi-conducting polymers are 

the basis for organic light-emitting diodes (LED's) and transistors [5] ; 

oligomers can be used in photovoltaic devices [33] . 

An important open question about these chemical structures is 
whether and how it is possible to extrapolate the electronic struc- 
ture of oligomers with increasing number m of segments to that of 
the polymer. The asymptotics of the band gap AE{m) of oligomers 
with length m increasing to infinity is subject to ongoing debate [4], [T], 
I17L I27L I28L I32L 143] . It is known both from experiments and quantum 
chemical calculations that the convergence behavior of the band gap 
of oligomers with increasing length to the band gap of the polymer 
often deviates from AE{m) — AE{oo) oc 1/m. One might expect the 
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(a) (b) 

Benzene Naphthalene 




(c) /^-carotene 

Figure 1. Molecules with conjugated vr-systems. Car- 
bon atoms are conventionally points connected by bonds. 
Hydrogen atoms are not displayed. 

latter relation based on the behavior of a linear chain of length m in 
the limit m — > cxd. The latter can be calculated straightforwardly and 
has been observed experimentally [2] (the finite chain corresponds to 
a C2mH2m+2 poly-acetyleue oligomer, PAm). 

We will use the following notation for oligomer /polymer molecules: 
the polymer is abbreviated by 2-5 letters related to its full name, e.g. 
PA (poly-acetylene), PPP (poly-para-phenylene). A number m follow- 
ing this abbreviation refers to an oligomer with m copies of the corre- 
sponding monomer, e.g. PPP5 is the para-phenylene oligomer C30H30. 
It consists of 5 connected copies of the CeHg monomer benzene. 

In this paper we discuss the gap AE{m) and its convergence in 
the polymer limit (m — > oo)for quantum chemical models that reduce 
to solving an eigenvalue problem on either combinatorial or metric 
graphs. Some background on these models will be given in Section [2l 
In Section Owe will define the spectral problems on graphs and describe 
how they are related to oligomers. We will also explicitly solve the 
simplest polymer (poly-acetylene in either ring or chain formation) as 
an example in Section 13.11 and summarize the (well-known) Floquet- 
Bloch theory for the infinite polymer in Section [321 Eventually we will 
present and discuss numerical results for a number of relevant polymers 
in Section H] and compare those to experimental data. 

2. Conjugated double bonds: beyond valence bond theory 

Conventionally, organic molecules are described by covalent bonds. 
According to the valence-bond (VB) theory the outer shell (valence) 
electrons spatially rearrange themselves such that a pair of electrons 
is shared by two atoms while inner shell electrons remain assigned to 
individual atoms of the molecule. A single covalent bond consists of one 
pair of electrons. Double and triple bonds between two atoms consist 
correspondingly of two or three pairs of electrons shared between the 
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(c) PPf 

Figure 2. Conjugated polymers poly- acetylene, poly- 
par a-phenylene, poly-pentafulvene. 

atoms. Recall that one single carbon atom has four valence electrons 
and it forms hence four bonds. Some of these may constitute a double 
or triple bond. A sequence of alternating single- and double C-C bonds 
is called conjugated, see Fig. [T] for examples. 

The VB assumption fails to explain the properties of compounds 
with conjugated sequences already on a qualitative level. Absorption 
spectra, luminescence and conductivity indicate that the spacing be- 
tween energy levels in such molecules is much lower than in similar 
non-conjugated counterparts. This strongly indicates that some elec- 
trons in these molecules are delocalized over a much larger region than 
one single or double bond, in contrast to the VB theory. 

The more general Molecular Orbital (MO) Theory overcomes the 
shortcomings of the VB theory |40] . It postulates that valence elec- 
trons recombine and form molecular orbitals to minimize the total en- 
ergy of the entire molecule. The following two models have been put 
forward to describe conjugated systems in the framework of MO The- 
ory. Both neglect the repulsion and correlation of electrons. Erich 
Hiickel suggested the following scheme for molecules with alternating 
double bonds [T3l [T4l fT5l [T6] : 

i. Single covalent bonds between adjacent atoms are formed and 
they determine the geometry of the molecule. 

ii. The remaining electrons (also called n- electrons, in a conju- 
gated system there is one 7r-electron per carbon atom) and 
their wave functions are subject to energy minimization. This 
reduces to an eigenvalue problem on a combinatorial graph af- 
ter some further simplifying assumptions. This spectral prob- 
lem is referred to as the Hiickel Molecular Orbital (HMO) The- 
ory. 

The alternative to the HMO approach is based on the assumption that 
the electrons of a conjugated system are strongly bound to the pos- 
itively charged carbon-backbone [38] . which is modeled as a graph 
I26L 137] . Interestingly, in the original formulation Hans Kuhn already 
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considered atoms as scattering centers for the electron wave functions. 
This is the Free Electron Model (FE). This model can be reduced (after 
separation of variables) to an eigenvalue problem on a metric graph for 
the Schrodinger operator without a potential (hence the 'free electron' 
in the name), u The assumptions of this model are only plausible for 
large systems. 

For finite molecules both, the MO and the FE model, yield dis- 
crete energy levels En with corresponding eigenfunctions (molecular 
orbitals). Owing to spin, each of these levels can be occupied by two 
electrons (applying Hund's rule and the Pauli exclusion principle). In a 
system with M alternating double bonds there are 2M vr-electrons. In 
the ground state of the system, the M lowest molecular orbitals are oc- 
cupied. The energy of the highest occupied molecular orbital {HOMO) 
is thus Ehomo = Em- The energy of the lowest unoccupied molecular 
orbital (LUMO) is Elumo = -E'm+i- 

The difference AE = Elumo — Ehomo is the smallest energy package 
a molecule can absorb. It is directly related to the color of a mate- 
rial (i.e. to the highest absorbed wave length) and to the breakdown 
voltage of semi-conducting materials (since this is proportional to the 
band-gap). In Ref. |42] band gaps of several semi-conductive polymers 
have been studied qualitatively using HMO theory. 

The degeneracy of Em can determine chemical properties of a com- 
pound. If Em is a multiple, completely occupied energy level, then 
the compound is chemically stable and inert. The classical example is 
nitrogen. Conversely, if Em is multiple and not completely occupied 
[SOMO, where "S" stands for "singly"), then the compound is chemi- 
cally reactive. The classical example is oxygen. The existence of such 
orbitals in conjugated hydrocarbons is related to the question whether 
is an eigenvalue of the adjacency matrix of the underlying graph, see 
[HI Chapter 8] which is dedicated to this topic. 

One advantage of the HMO theory is that absolute energy levels 
can be calculated. As a consequence it is possible to investigate the 
stability of compounds. Moreover, the underlying MO theory is not 
restricted to conjugated systems. 

On the other side, the FE model has the advantage that no parame- 
ter fitting is involved. Moreover, wave functions and electron densities 
can easier be calculated in the FE model. The latter information can 
then be used to infer statements about bond lengths and effect of het- 
eroatoms 



Arguments of [37j are incorrect for MO's localized at an atom (eigenfunctions 
of the 3D-Schr6dinger operator localized at a vertex) |21] . Such states would cor- 
respond according to the MO Theory to unbound electron pairs and are, therefore, 
irrelevant for realistic hydrocarbons. 
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3. Mathematical description of polymers by combinatorial 
and metric graph models 

For a given molecule with a conjugated fragment and Nc carbon 
atoms, it is straightforward to define the underlying graph G. The 
vertices Vi (where 2 = 1,2,..., Nc) of the graph are the carbon atoms 
and the undirected bonds h of the graph are the chemical bonds be- 
tween them. For the undirected bond connecting the vertices Vi and 
Vj we write hij = bji. Note that on the graph G all bonds are simple, 
disregarding their chemical multiplicity. We write i ~ j if the vertices 
(atoms) Vi and Vj are connected. The undirected (metric) graph is 
fully characterized by its Nc x Nc symmetrical bond matrix B, where 
Bij = Bji is the bond length between the two carbon atoms (vertices) 
Vi and Vj if they are connected, and otherwise. In the following all 
bonds in a purely hydrocarbon system are assumed identical, unless 
the opposite is stated. By replacing all bond lengths in B with unity 
one obtains the Nc x Nc connectivity matrix C. 

HMO theory in its original formulation solves the eigenvalue prob- 
lem 

(3.1) Cu = Xu. 

Since C is a real symmetric matrix, it has Nc real eigenvalues Ai, . . . , X^c- 
The electronic energy levels En in a molecule are related to the eigen- 
values Xn of the connectivity matrix by 

(3.2) En = Xnf3-a, n = l,...,Nc 

where /S and a are constant parameters with units of energy. The value 
of a is not relevant for pure hydrocarbon compounds since it is just 
a constant shift. The parameter (3 is defined as a resonance integral, 
but is usually fitted to experimental data. In order to keep results 
consistent with the FE model, we choose /? = 3.05eV as in Ref. |41] . 
It is important to note that the parameters a and (3 are not potentials 
in the physical sense, and that the eigenvectors Un are not the values 
of the wave function at the vertices. 

To formulate the FE model, we need to define a real coordinate 
< Xij < Bij along each undirected bond bij. For i < j we set 
Xij = at one end of the bond where it is connected to the vertex Vi 
and Xij = Bij at the other end connected to Vj. We will also use the 
notation Xij = Vi for the value of Xij at vertex Vi. The metric graph 
is quantized by defining a self-adjoint Schrodinger operator for a wave 
function $(x) on the metric graph. 

This implies that the wave function <l>(x) must satisfy the one-electron 
free Schrodinger equation on each bond bij 

(3.3) __^-^»(,,^) = s*(,.^) 
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with Plank's constant h, electron mass m^ and energy E. 
To obtain a self-adjoint Schrodinger operator one has to add certain 
conditions on the wave functions at the vertices [H], I18L 119] . We will 
impose the Neumann conditions (also known as Kirchhoff conditions) 
that the wave function is continuous at each vertex 

(3.4) for all j ~ i lim ^(xij) =$(fj) 

and, additionally that the sum over all outgoing derivatives of the wave- 
functions on all bonds connected to a vertex vanishes 

(3.5) y: ^5d^*(^^^ 

— LJ ~ t. 

where w* = 1 if i < j and w* = — 1 otherwise. 

Solving the Schrodinger Equation 03. 3p with vertex conditions 03.41 - 
13. 5p we obtain discrete energy levels En and wave functions (MOs) 

Additional potentials on the bonds have been introduced in Eq. 03.31) 
to improve the agreement with experimental values |26[ , I37j . One 
may also add vertex potentialqj to get a more realistic description (see 
[9l ll9] ). Any such generalization adds some fitting parameters. In this 
work we are mainly interested in more fundamental questions concern- 
ing the convergence behavior for general polymers and we thus stick to 
the simplest model and no fitting parameters. 

We scale the bond lengths of the graph to the average C-C bond 
length of Lc-c = 1.4A and -^ to 1. The eigenvalues /!„ of the scaled 
problem correspond to electronic energy levels as E^ = fJ^n^ with e = 
T^^ — = 1.95eV. 

In order to find the eigenvalues numerically, one has to solve the 
secular equation (see [9], 119] for details). However, throughout this 
work we assume equal bond lengths. Then the following identity is 
valid |34j . Define the diagonal valency matrix V (i.e. Vi^i is the number 
of bonds at the vertex Vi). If /i is an eigenvalue of the (scaled) FE 
problem, then fl = cos ^/JI is an eigenvalue with the same multiplicity 
of the generalized eigenvalue problem 

(3.6) Cu = flVu. 

We adopt the following notation for oligomers and polymers, see 
also Fig. [3l A monomer X = XI is described by a graph that we will 
call Gi. The number of carbon atoms in the monomer will be denoted 
by Nc and is equal to the number of vertices in Gi, and N will be used 
for the number of double bonds in the monomer. Usually, Nc = 2N. 



A vertex potential is equivalent to changing the inner condition (j3.5p at a 
vertex while keeping the continuity condition p.4p 
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Figure 3. Construction of the oligomer PPfm, chemical 
notation is used to represent the graph. The fifth carbon 
atom of the first monomer in the chain is connected to the 
third atom of the second one etc. In the graph notation 
this means v^ ~ w|; the index fe is 5 and the index e is 3. 



The monomer has two atoms at which bonds to other monomers of 
the same type are created. Let these atoms be represented by the 
vertices Vb,Ve, so that two consecutive monomers are connected by a 
bond between the atom b of the first and e of the second monomer. 
Let us now consider m copies Gi, G2, G3, . . . , Gm of the monomer Gi. 
Denote the vertex Vi on Ga by vf. We can now construct the graph 
G"^ that represents the oligomer Xm as a union G"^ = IJ^i ^a with 
extra bonds between v^ and v^~^^, a = 1, . . . ,m — l. Obviously one has 
Gi = G\ The band gap of G"" is 

(3.7) AE{m) = Em+i - Em 

where M = Nm is the number of double bonds on the oligomer. 

The above construction assumes that adjacent monomers are con- 
nected via a single bond between two atoms. We only consider poly- 
mers of this type here. Note that the construction can easily be gener- 
alized to incorporate many bonds between adjacent monomers. 

3.1. Examples: rings and chains. Consider a ring C2mH2m con- 
sisting oi M = m alternating double bonds. It has Nc = 2 carbon 
atoms per "monomer", each with A^ = 1 (trivially) alternating double 
bond. The HMO-eigenvalues are 



(3. 



,mr , 



A„ = 2cos( — ), n= 1,2, ...,2A^ 
m 



while the FE eigenvalues are 



(3.9) 



/in 



n^TT^ 



m^ 



n 



0,±1,±2, 



We see that both models predict qualitatively the same energy scheme. 
The ground state is simple, higher states are double (except for A2m, 
which is not important). Thus, Em and Em+i are both double. For odd 
m, Ehomo is completely filled, yielding an increased chemical stabihty, 
while for even m it is not complete leading to an increased reactivity. 
This is the famous Hiickel rule. 
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Consider now a chain C2mH2m+2 with M = m alternating double 
bonds (polyenes). The limit ?ti — >■ cx) is the most trivial example for a 
conjugated polymer and is called poly-acetylene (PA). It was the first 
polymer for which electric conductivity was shown and studied. 

The HMO and FE eigenvalues are 

(3.10) ^" = ^''°H2^^^TTJ' ^ = l'2'---'2m 

71. TT 

(3.11) I^n = j7^ rr^, n= 0,1,2,... 

respectively u. Now all energy levels a single. The HOMO-LUMO 
difference in HMO is 

(3.12) AEHMoim) = E^+i - Em = -4/3 sin ( -^^ 

\Am + 2 

and in FE 

/ 2m + 1 

(3.13) AEpEim) = Em+i - Em = en' i——- 

We see that both models predict AE{oo) = and AE[m) — 
AE{oo) oc — for large m. In both models, the discrete electronic 
levels become one single band in the m — > oo limit. In practice, the 
band gap of PA is significant (1.8eV), but the convergence behavior 
AE{m) — AE{oo) is almost linear in 1/m, with the slope close to the 
FE and HMO case, see Section HI 

3.2. Polymer limit. In the limit m — > oo the underlying graph 
G°° is periodic and Floquet-Bloch Theory can be applied [T], IllL I20L 
[m [Ml [211 [HE]. Consider the HMO eigenvalue problem (the FE case 
is analogous and will not be discussed in detail) 

(3.14) Cgoo^ = A*. 

All physically relevant eigenvector candidates (bound states) can be 
found as a product of a periodic function and a phase shift: 

(3.15) <i,(vr)=e^>^Mj(vi), 

where we introduced the quasi-momentum k G [— 7r,7r]. For a multiple 
eigenvalue A, Floquet-Bloch eigenf unctions span the eigenspace. The 
eigenvalue problem on G°° reduces to the problem on a single monomer 

(3.16) (C + F(A;))^ = A(A;)^ 

where F[k) is a Nc x Nc matrix with only two non-zero entries at 
the connection between monomers Feb{k) = e^^ and F},e{k) = e~^^. 



In the case of the free electron model (metric graphs) it seems rather unnatural 
to assume that the MOs in a linear chain are limited by the end atoms. We thus 
add dangling bonds at vl and u™, so that the total length of the MO becomes 
2iV + 1 [m [H], yielding A„ ^ j£^ and AEMm) - 2^- 
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Since C + F{k) is a Hermitian matrix, it has A''^; real eigenvalues 
Ai(/c), . . . , AAr^(/c). Each of the branches \r{k) is a continuous func- 
tion of k. The set IJo<fc<7r •^»'(^) ^^ ^^^ ^^^ band. The iVth band is 
called the valence hand while the (A^ + l)th band is the conduction 
hand. 

The spectrum of Cg°° is the union of all bands 

(3.17) a(CGoc)= U U \rik). 

0<k<T7r=l,...,2N 

Due to the one-dimensional periodicity and the fact that the monomers 
are connected by one single bond, the extrema of Xr{k) are always at 
k = or k = IT. Consider the characteristic polynomial ^ of C -|- F[k) 
as a function of k and A 

(3.18) ^{k,X) = det{C + F{k)-XlNc) 

where Inc is the unit matrix. Eigenvalues are given by C,{k, A) = 0, 
and hence the Implicit Function Theorem can be applied: if g| 7^ in 
a neighborhood of k, A with C,{k, A) = 0, then A(fc) is a well-defined, 
differentiable function in this neighborhood with 

dX fd^\ fdC 



^^•^^^ dk \dXr~\dk 

One easily finds that 

(3.20) 1^ = te'^i-iy-'M^^t, - te-^^i-iy^'Mt^e = 
ok 

~2{-lY+''lm{e'''M,^b) = Im [e''' (mo(A) + mi(A)e-*'=)] = mo(A)sinA; 

where M^^b is the e, 6-minor of C + F{k) — XI^^, i.e. the determinant 
of the matrix after deleting the eth row and 6th column; rriQ, mi are 
some polynomials in A. The first equation follows from the rule for 
determinant derivatives; the second one from the fact that M^^b = Mb^e] 
the third one because Me^ has only one entry which depends on k; and 
the last one because mi (A) is real. We obtain therefore 

(3.21) ak, X) = po(A) cos(fc) + pi(A), II = p[,(A) cos(A;) + ^(A) 

with some polynomials po (= —mo) and pi. Obviously, X{k) = X{—k) 
and A determines fc G [0,7r] uniquely iff Po does not vanish. Now three 
cases are possible for an eigenvalue A: 

(1) po{X) = pi(A) = 0. Then C,{k, A) does not depend on k, and A 
is a degenerate band A(/c) = A. And, conversely, if A belongs 
to a non-degenerate band, then po{X) 7^ 0. 

(2) g|(A) = 0. This means that A is a multiple eigenvalue, i.e. 
belongs to two (or more) bands. In this case, the two bands 
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join. Assume namely, that two non-degenerate bands over- 
lap, i.e. there exists an interval (A_, A+) of "double" eigenval- 
ues. Then, according to Eq. fl3.2ip . almost everywhere on this 
interval would hold cos{k) = Pi{\)/po{\) and p'q{\)pi{\) + 
p'iWpoW = 0. Therefore, both po and pi must be constant, 
which is only possible for empty graphs. 
(3) Otherwise, po{\) 7^ on the entire band, and the extrema of 
X{k) within this band are given by the zeros of Eq. fl3.20p . i.e. 
integer multiples of vr. 

Thus, in all cases the band edges are contained in the eigenvalues of 
C + F(0) and C + F{7i) {e"-^ = ±1). Note that in the second case no 
statement is made for degenerate bands. Remarkably, if F{k) has more 
than two non-zero entries, i.e. the monomers are connected by more 
than one bond, then the arguments of Eg. (13.201) become incorrect, and 
a counterexample to the statement about band edges can be found in 

m- 

To study and compare spectral bands of oligomers and polymers, 
it is useful to define for a real number x the spectral counting function 
as the number of eigenvalues per monomer, which are smaller than x. 
For an oligomer this is given by 



(3.22) 



cr [X 



{n|A„_i <x < \n} 



m 



where the eigenvalues Aq, Ai, . . . are put in increasing order and A_i = 
—00. The spectral counting function of the polymer is obtained by 
taking the limit ?Ti — * 00. 

Intervals where a increases continuously correspond to bands; inter- 
vals where it is constant correspond to band gaps. Jumps are de- 
generated bands: e.g. if Gi has an eigenvalue A with an eigenstate 
which vanishes at the vertices Vh and v^ then the oligomer will have 
an ?Ti-fold degenerate eigenvalue with quasi-momentum fc = and 
Ai(0) = A2(0) = ■ ■ ■ = Am(0) = A with eigenstates which are local- 
ized on one monomer. Such eigenstates do not satisfy Eq. fl3.15p . but 
lie in the space spanned by Floquet-Bloch eigenfunctions. 
Consider PPf, see Fig. [3l The matrices C and F{k) are 

(3.23) 



C 



/ 1 1 \ 




/ 





\ 


10 10 













10 10 
10 11 


; F{k) = 










e'^ 



10 10 







g-ifc 





\ 1 / 




i^ 





y 



The calculated (t'^^q 



is shown in Fig. H] together with the band edges 
from Floquet-Bloch Theory. The first two bands are [—2.34, —2] and 
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PPf HMO 




-1,56-1,41 



-0.47 (dble) 



1 (dble) 1,41 1,81 



Figure 4. agJJJo and Floquet-Bloch Theory of PPf. The 
twelve band edges arranged in increasing order —2.34, 
-2, -1.56, -1.41, -0.47, 0, 0, 1, 1, 1.41, 1.81, 2.56 are 
displayed on the horizontal axis. 

[—1.56, —1.41]. Since and 1 appear twice, there are no gaps between 
the third, fourth and fifth bands. The last band is then [1.81,2.56]. 
Knowing that is a double eigenvalue of C + F{0) while 1 is a simple 
eigenvalue of C + F(0) and C + F{7i), we conclude that there is a jump 
at 1 and no jump at 0. 



4. Examples and comparison with experiment 

It is important to note that experimentally determined band gaps 
systematically differ from the ones predicted by simple theories. The 
differences are due to (i) intermolecular phenomena such as electronic 
repulsion, electron correlation, spatial geometry, non-conjugated side 
chains and (ii) intermolecular effects such as solvent or solid-state shifts 
and impurities. These effects are approximately constant (~ leV) for 
a given sequence of oligomers at the same experimental conditions 
[7]. The band gap of a (real) material is influenced by the environ- 
ment, therefore different measurements on the same compound may 
give slightly different values depending on experimental conditions. 

We consider the following polymers: PA [Nc = 2, N = 1), PPf 
(poly-pentafulvene, A^^ = 6, iV = 3), PPP {Nc = 6, N = 3), PMP 
( poly-met a-phenylene, Nq = Q, N = 3), PPV (poly-phenylene-vinylene, 
Nc = 8, iV = 4) and PmPV (poly-meta-phenylene-vinylene, Nc = 8, 
N = 4), see also Fig. [2l The first one is the alternating chain and 
has already been discussed in Section 13.11 The second one, PPf rep- 
resents the large number of 5-ring based polymers. To the best of our 
knowledge, it is subject to theoretical studies |12U35] as a possible low- 
band-gap polymer, but has not been synthesized yet. The other four 
are based on benzene (6-ring) as monomer but connected differently. 
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Figure 5. Band gaps of oligomers estimated by HMO 
and FE models compared with experimental values. On 
the horizontal axis is the inverse number of monomers. 
HMO-values are represented by circles; FE-values by 
stars. Sources: all TD DFT calculations [29] : PA [2]; 
PPP [30], [31]; PPV [6], [8], [7]; PmPV [ID]. 



PMP has a 1 — 3 connection [meta) whereas PPP has a 1 — 4 con- 
nection {para). PPV and PmPV have additionally two carbon atoms 
between the rings, i.e. the connection length is 3 and A^ = 4. Note 
that PPP2=PMP2 and PPV2=PmPV2. 

As we can see from Fig. 0, both models predict similar trends and 
values. The estimated band gaps are mostly within leV below the 
experimental values. PPf is the only polymer showing a PA-like be- 
havior: the band gap converging to as 1/m. This can be explained 
by symmetry arguments [42] . The non-linear (in 1/m) convergence of 
other oligomer series is in line with experimental observations. The 
band gap decreases by 1.5 — 2eV between m = 2 and ?7i = oo in se- 
ries of para-connected oligomers PPPtti and PPVrra, whereas for the 
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PAFE 



PPf FE 




Figure 6. Spectral counting function a^^ of PA, PPf. 
PPP, PMP, PPV and PmPV in the FE model 



meta-connected PMPm and PmPVm this decrease is only ^ 0.5eV, 
in a good agreement with available experimental data. We have also 
compared our results with Time Dependent Density Functional Theory 
(TD DFT) calculation [29], data represented by triangles. These data 
compare better with measurements than the FE or HMO predictions, 
but provide less information about trends. 

Another crucial property for conductivity is the width of the valence 
and conduction bands. We have calculated cthmo ^^^ "^fe ^^^ ^^^ ^^^ 
above polymers {m = 200) in Figs. [6] and [71 

As discussed above, the band gaps of PA and PPf vanish. We can 
see that the para-connected polymers have much wider valence and 
conduction bands than their meta-connected counterparts. Together 
with the smaller band gaps this explains why PPP and PPV are good 
conductors while PMP and PmPV are not 
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FlGURE 7. Spectral counting function (Thmo o^ P^' PP^' 
PPP, PMP, PPV and PmPV in the HMO model 



5. Conclusions 

We have studied conducting polymers from the point of view of 
Hiickel Molecular Orbital Theory and Free Electron Model, the two 
models directly related to the spectral graph theory. Qualitative pre- 
dictions match the experimental reality. Floquet-Bloch theory for pe- 
riodic graphs provides means for an easy study of band properties of 
polymers. 

Quantitative comparison shows that actual band gaps can be rea- 
sonably estimated by either of the models when correcting for in- 
tramolecular and intermolecular shifts. The absence of any parameter 
fitting makes this result particularly exciting. The convergence behav- 
ior AE{m) — AE{oo) of band gaps of oligomer series is well reproduced. 
Therefore, it is possible to extract the band gap of the polymer from the 
data for finite values of m with sufficient accuracy using simple graph 
models. It would be promising for practical applications to derive an 
analytic expression for AE{m). 
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